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Abstract 

We present a simple numerical procedure for calculating the irrotational hydrodynamic flow in 
a helium solvation structure around a spherical solute in linear motion through superfluid helium. 
The calculation requires only the radial helium density around the impurity as input. From the 
resulting irrotational flow, the helium contribution to the effective mass of the solute is calculated. 
The results for alkali cations are compared to recent many-body Variational Monte Carlo (VMC) 
calculations by M. Buzzacchi, D. E. Galli, and L. Reatto (Phys. Rev. B., 64 094512 (2001)). The 
helium contribution to the effective masses calculated by the two methods are 12.9(4.6) versus 9.4 
u for Li+, 48.2(5.6) versus 52.1 u for Na+, 69.6(4.8) versus 70.1 u for K+, and 6.4(8.8) versus 6.8 for 
Cs + , with the VMC result listed first (with one a statistical error estimate) and the hydrodynamic 
result listed second. For the cases of Na + and K + , the hydrodynamic calculation treated the 
first helium solvation shell as a rigid solid, as suggested by the VMC calculations; treating the 
first solvent layer as part of the superfluid lead to considerable underestimate, ~ 50%, of the mass 
increase in both cases. In all cases, the agreement of the two results are in within the error estimate 
of the VMC calculation, demonstrating the accuracy of hydrodynamic treatment of helium motion 
even on the atomic length scale. 
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The study of atomic ions as microscopic probes of superfluidity in helium-4 has a rich 
history Recent progress on both the injection of neutral metal atoms into bulk liquid 
helium [fj, |3|] and the study of doped helium nanodroplets ||] have caused a renewed in- 
terest in the dynamical behavior of solutes in this unique quantum liquid. In particular, 
concerning the rotational degrees of freedom, here have recently been several papers that 
try to provide a microscopic explaination for the observed sizable fractional increase in the 
rotational moments of inertia of all but the fastest rotors when dissolved in superfluid he- 
lium I, | | I, 0. In all of the proposed models, the extra moment of inertia arises 
from helium kinetic energy induced by rotation of the molecule, but there is disagreement 
about the way this should be calculated, and about the physical description of the helium 
motion [p"Tf . 

A closely related problem is that of the change in the effective translational mass of 
a solute in superfluid helium, which also arise from helium kinetic energy induced by the 
requirement that the helium solvation structure must 'follow' a moving impurity. In the 



quantum hydrodynamic model |q, 12], one calculates the classical velocity potential that 



describes the irrotational flow which maintains a constant helium solvation structure in the 



frame of a moving impurity. By a theorem of Lord Kelvin || |r5| , this irrotational flow will 
provide the minimum kinetic energy flow that satisfies the equation of continuity. A simple 
example is the case of a moving spherical "hole" in the liquid, for which the hydrodynamic 
flow is a dipole field that carries a kinetic energy equal to one half of that of the liquid mass 
displaced by the sphere moving at the velocity of the hole ||13|| . For the general case of an 
inhomogeneous density, Barrera and Baym [14| have presented a solution to the equation 
of continuity, based upon a transformation of dipole solution. However, this transformed 
solution is not irrotational, and if fact has a continuously varying vorticity, and thus is not 
appropriate for flow in a superfluid, where vorticity must be quantized. 

This paper presents a general numerical scheme for finding the irrotational solution of 
the equation of continuity given a solvation density around a moving spherical solute. The 
hydrodynamic equation can be solved by separation of variables, and reduces to one di- 
mensional quadrature of the radial homogeneous and inhomogeneous equations, allowing 
the mixed boundary conditions to be satisfied without need for iteration. In the limit of 
infinitesimal solute velocity, the given solution can be shown [|EJ] to provide a variationally 
optimized helium ground state wave function, assuming a one-body phase function. As such, 
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it should provide a rigorous lower bound on the true increase in effective mass. 

In order to test the quantitative accuracy of this approximation, we have compared the 
calculated mass increase for Li + , Na + , K + , and Cs + with those calculated by Buzzacchi, 
Galli, and Reatto |T5j using a Variational Monte Carlo (VMC) treatment of the explicit 



many-body problem [ 03 , . The input to the hydrodynamic calculations were the helium 
radial densities calculated by these same authors using the same method; thus the compar- 
ison provides a direct test of the accuracy of the one-particle hydrodynamic treatment on 
atomic length scales. 

I. CALCULATION METHOD 

Consider an atomic solute in superfluid He that has a solvation structure with radial 
number density p(r). We assume that the solute is moving with velocity V (in a direction 
we take as the z axis), and that the solvation structure adiabatically moves with the atom. 
This generates an irrotational flow in the helium of velocity v = — V0, where is the 
velocity potential, which must satisfy the equation of continuity: 

V(pV0) = f t = -(Vp) • (Vz) (1) 

If we write <fi = <p(r) V cos(9) where r is the distance from the impurity atom and 9 is the 
angle with the z axis, we find that the hydrodynamic equation of continuity is solved if 
satisfies the following radial equation: 

^ + lf.^ + (^i)(f)^(^i) (2) 

^ r z r d r r * \ ar J \ar I \ ar J 

At long range from the atom, the density must approach the bulk value, p e , and — > A/r 2 . 
On the inner wall of the solvation structure, r ; (where the helium density vanishes rapidly), 

— ► — 1. The general solution to the inhomogeneous equation ||] can be written as a sum of 
any inhomogeneous solution plus any linear combination of the homogeneous solutions. The 
homogeneous equation is found by setting the right hand side of Eq. |2| equal to zero. The 
inhomogeneous solution that satisfies the boundary conditions was found as follows. At large 
r D , inhomogeneous and homogeneous solutions were started assuming the asymptotic form. 
The value of A = rf/2, which is the correct value for the uniform density case, was used 
to start the solutions. The two solutions are numerically integrated until r = r ; , at which 



point the derivatives of the homogeneous solution, </>h and the inhomogeneous solution, (j)^ 
are used to determine the constant B, equal to: 



1 + 



B 



dr 



(3) 



The inhomogeneous solution satisfying the boundary conditions is: 



<t>( r ) = 4>mh(r) +B(f> h (r) 



(4) 



The hydrodynamic kinetic energy is found by integrating the helium kinetic energy, which 
is proportional to V 2 . This allows us to define a hydrodynamic mass, M h , by 

M h 



Mi 



He 



V~ 2 J p\V<p\ 2 dV 
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p(r) 



dr 



+ 



20(r) 



r dr 



(5) 



This integral is evaluated from the numerical solution over the domain r ; < r < r Q . If we 
assume that the solution for r > r G is given by the asymptotic form, then this gives an 
additional contribution of (8ir /3) Mn c p e <f)(r ) 2 r to the integral defining M^. It has been 
checked that this definition of gives the correct value of one half the displaced helium 
mass for the case of a hole in helium of uniform density. 

It is also possible to use the hydrodynamic equation to transform the integral for the 
effective mass to give: 
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He 
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dp 
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r ° ~ f dp 
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(6) 



There is no contribution to the volume integral in the region r > r D because by assumption, 
(clr) = i n this region. The two estimates for Mn e need agree only if is a solution of the 
continuum hydrodynamic equation, and thus a comparison between them provides a test of 
the convergence of the numerical solution and the size of spacing used for integration. 
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II. APPLICATION TO ALKALI CATIONS 



Figure 1 shows the helium solvation densities around Li + , Na + , K + , and Cs + , as cal- 
culated by VMC [15| using a trial function of the "shadow function" form. This tech- 
nique treats solids, liquids, and solid-liquid mixtures with a single functional form citePed- 
eriva94,Duminuco00. It can be noted that the solvation density goes almost to zero between 
the first and second solvent layers for K + , and is highly structured in the case of Na + . The 
VMC calculations find that the first solvation layer has a highly solid-like order for the Na + 
and K + cases with little to no exchange of these atoms between solvent layers [|I5|, [T7| . This 
suggests that the first solvation layer not be treated as part of the fluid but as a fixed mass 



that moves rigidly with the cation, as in the snowball model of Atkins ljl8|| . In contrast, the 
VMC find substantial mobility between solvation layers for Li + and Cs + , suggesting that in 
these cases that even the highly compacted first solvation layer should be treated as part of 
the superfluid, and thus be treated as part of the hydrodynamic flow. 

Table | contains comparisons of the effective mass for each cation calculated using a 
hydrodynamic treatment and estimated from the VMC calculations. It is seen that in all 
four cases, the agreement of the two estimates is excellent, being within the VMC statistical 
error estimate. For the cases of Na + and K + , the first solvation shell, with 10 and 12 
helium atoms respectively, was treated as a rigid solid and the hydrodynamic calculations 
were begun at the minimum density point between the first and second solvent shells. In 
both these cases, treating the entire density with the hydrodynamic approach yielded a 
substantial underestimate for the translational mass; 21.7 u for Na + and 36.0 u for K + . The 
hydrodynamic calculations are vastly less computationally expensive than the many-body 
treatment [ |H| . 

It is useful to compare our present results with that of the widely used model of the 



cation "snowball" due to Atkins |18fl . In this electrostrictive model, the helium is treated 
as a continuum dielectric material, whose density is increased near the cation due to the 
ion-induced dipole interaction. For a radius less than b, on the order of 5-6 A (whose value 
depends upon the helium liquid-solid surface tension assumed but is independent of the 
specific singly charged cation |J), the helium is predicted to form a solid "snowball" that 
moves rigidly with the ion. This snowball contributes a mass of ~ 150 u to the ion effective 
mass n. In addition, there is a hydrodynamic contribution to the mass, expected to be 
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on the order of the hard sphere value mns = ^-b 3 p e mn e ~ 40 u. Taking into account the 
increased helium density for r slightly larger than b, using the model used by Barrera and 
Baym [14]] (p — p e = p e X(b/r) 4 with A = 0.186), the above hydrodynamic treatment predicts 
the hydrodynamic contribution to the mass to be 0.932mHs, which can be compared to 
the value 0.97mns reported by Barrera and Baym |T^] for their proposed velocity solution 
which is not irrotational. In agreement with the Kelvin minimum energy principle [^] , their 
solution is higher in kinetic energy and thus predicts a higher mass. Comparision with the 
both the VMC and hydrodynamic results show that the effective masses, even for the case 
of rigid first solvation shells, is considerably less than those predicted by the snowball model 
of Atkins. 



III. CONCLUSIONS 

The present work demonstrates that the hydrodynamic treatment of the linear motion of 
a solute through superfluid helium predicts the solvent contributions to the effective mass of 
the impurity in quantitative agreement with more exact many-body approaches, yet requires 
only a trivial additional computational cost once the solvent density has been calculated. 
This applies to helium density well inside the predicted 5-6 A radius of the liquid-solid 
surface in the "snowball" model of Atkins Hl8~| . However, the highly ordered first solvent 
layer around some ions must be treated as a solid that rigidly moves with the ion. The 
present results compliment our recent hydrodynamic calculations of solvent contributions 
to the moments of inertia of molecules solvated in helium, which were found to be in good 
agreement with experiment ||. The present results have been directly compared to higher 
levels of theory and thus provide a more critical test of the hydrodynamic model, since 
uncertainties in the solute-helium potentials do not enter. We therefore demonstrate that 
the hydrodynamic treatment of superfluid helium motion can be accurate on the atomic 
scale, directly refuting recent criticisms of its use 0. 
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(this work). For the cases of Na + and K + , the hydrodynamic calculation treated the first helium 
solvation shell as crystalline, as predicted by the VMC calculations. 



Cation VMC Hydrodynamic 

Li+ 12.9(4.6) 9.4 

Na+ 48.2(5.6) 52.1 

K+ 69.6(4.8) 70.1 

Cs+ 6.4(8.8) 6.8 
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